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I. INTRODUCTION 

Motivated by the discovery that quantum algorithms 
can provide exponential gains for solving certain com- 
putational problems, numerous proposals have been 
advanced for experimental realization of a quantum 
computer While a useful processor remains far 

off, ground breaking experiments have demonstrated 
controlled evolution of a few qubits and implemented 
basic quantum algorithms for computation and error 
correction^, 0, IE IE 0, 0. Among the range of phys- 
ical systems that have been identified as candidates for 
implementing quantum logic the most extensive labora- 
tory results so far have been obtained with cold trapped 
ions £j and nuclear magnetic resonance (NMR) in macro- 
scopic samples [lfjl |llj. 

Within the last few years neutral atoms have emerged 
as a possible route to experimental quantum logic. 
The most obvious distinguishing feature between neu- 
tral atom and trapped ion schemes is the absence of 
strong Coulomb forces in the former. Coulomb forces 
between ions couple strongly the motional degrees of 
freedom. This can be utilized to entangle any two of 
a linear string of ions, as was first elucidated in the 
work of Cirac and Zoller^^ an d demonstrated experi- 
mentally in Boulder and Innsbruck 0|. The lack of 
a strong Coulomb interaction in neutral atoms is advan- 
tageous as regards decoherence, since coupling to stray 
fields is weaker for atoms than for ions. The draw- 
back, and indeed the central difficulty in constructing 
a large scale quantum processor, is the need for strong 
qubit to qubit coupling, while maintaining weak coupling 
to the environment. Neutral atom coupling based on 
ground state collisions, light mediated dipole-dipole cou- 
pling, and dipole-dipole coupling of highly excited Ry- 
dber g st ates have all been proposed in the last several 
yearsjH GJ, HE H HJ . In particular dipole-dipole 
coupling of Rydberg states provides a strong interaction 
suitable for the implementation of fast gates[2l|, and this 
paper is devoted to a detailed study of this approach. 

While the theoretical foundations of the Jaksch et al. 



Rydberg state dipole-dipole coupling approach to quan- 
tum logic have been presented[2l| the question of how to 
implement this scheme in a practical and scalable fash- 
ion has not been solved 22]. Regardless of how logical 
operations are to be performed, there are two primary 
obstacles that must be surmounted. The first is how to 
create a large number of trapping sites, and load a single 
atom into each site. This amounts to initialization of the 
quantum computer. The second difficulty is that in order 
to be useful for generic models of quantum computation 
the sites must be individually accessible for logical opera- 
tions, and state readout. In this paper we do not discuss 
the problem of creation and loading of a large number 
of addressable single atom sites. A number of possible 
solutions to these questions have been discussed in the 
literaturejH H HI |H |H |H |U Hj . 

Our goal here is to examine in detail the use of two 
closely spaced sites, each containing a neutral atom 
qubit, for high fidelity quantum operations. Far-off- 
resonance optical traps (FORTs) are defined by tightly 
focusing laser light in a set of chosen locations. Single 
atoms are loaded into the o ptic al traps after precooling 
in a magneto-optical trap |3lLl32l ]. Single qubit operations 
are performed using two-photon stimulated Raman tran- 
sitions, and a two-qubit conditional phase gate is realized 
using dipole-dipole coupling of atoms excited to high ly- 
ing Rydberg states[2lj. Qubit measurement is performed 
by counting resonance fluorescence photons. 

The ability to perform many operations with high fi- 
delity and low decoherence is a prerequisite for scaling up 
to a larger number of qubits. As will be shown in what 
follows our calculations lead to the conclusion that a set 
of logically complete qubit operations can be performed 
with high fidelity at MHz rates. This would suggest that 
qubit storage in optical traps with coherence times of 
much less than a second will be sufficient for large com- 
putations. However, the necessity of implementing error 
correction implies that a computation will also require a 
large number of state measurements, which are projected 
to be several orders of magnitude slower than the logical 
operations. We therefore examine closely the feasibility 
of Tl and T2 coherence times of several seconds in optical 
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traps. 

The remainder of this paper is structured as follows. 
In Sec. |n]we estimate the decoherence times for storage 
of individual atoms in FORTs. We specifically include 
the contributions due to collisions with hot background 
atoms (Sec. Ill Bl) . photon scattering from the trapping 
laser (Sec. Ill CJI . spin flips due to the trapping laser, and 
heating rates due to laser noise (Sec. I11D|) . Decoherence 
due to fluctuations in the trapping lasers is considered in 
Sec. Ill El and due to background fields in Sec. Ill Fl 

In Sec. IIIII we discuss the operation of single qubit 
gates based on two-photon stimulated Raman transi- 
tions. In particular we calculate decoherence proba- 
bilities due to excited state spontaneous emission (Sec. 
MI All and expected gate fidelities due to ac Stark 
shifts(Sec. MI Bp and motional effects (Sec. MI C|) . Leak- 
age out of the computational basis due to imperfect op- 
tical polarization is estimated in Sec. MI Dl and limita- 
tions imposed by the laser stability are estimated in Sec. 
MI El The ability to rapidly interrogate the atomic state 
is crucial to the usefulness of this approach. We discuss 
single atom state detection using collection of resonance 
fluorescence in Sec. MI Fl Included in Sec. MI Fl is a dis- 
cussion of heating during readout, and its amelioration 
using red-detuned molasses for the interrogation beams. 

In Sec. IIVI we discuss the implementation of a two- 
qubit conditional phase gate which can serve as a log- 
ical primitive for arbitrary computations. We consider 
two different regimes of operation: Rabi frequency large 
compared to dipole-dipole frequency shift (Sec. IIV A|) 
and dipole-dipole frequency shift large compared to Rabi 
frequency (Sec. IIV Bp . 
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TABLE I: Physical mechanisms limiting groundstate coher- 
ence of qubit basis states \a) — \F = l,m,F = 0} and 
|6> = \F = 2,m F = 0). Values listed are for P = 1 x 10" 10 
mbar, U m — 1 mK, T a — 50 /iK, A/ = 1.01 /im, and 
Wfo — 2.5 fj,m, and ato(uJf) = 114 A 3 . See text for details 
and additional parameters. 

II. OPTICAL TRAPS FOR SINGLE ATOMS 

In this section we recall some basic features of far 
off resonant optical traps. A number of distinct phys- 
ical mechanisms limit the coherence of atoms stored in 
FORT's. Some of these decohering effects are intrinsic to 
the operation of the FORT, and some are due to technical 
imperfections of the apparatus used. As shown in Table 
|U these mechanisms contribute to an effective decoher- 
ence time of diagonal (Tl) or off-diagonal (T2) density 
matrix elements. We discuss the physics behind each 
of these decohering mechanisms in the following subsec- 
tions. All numerical estimates of decoherence rates will 
be calculated for 87 Rb atoms using the parameters listed 
in Table HI 

A. FORT trap parameters 



Two-qubit operations involving Rydberg states have 
larger errors and higher decoherence rates than single 
qubit operations. We optimize the parameters of a phase 
gate in the two limits of weak dipole-dipole interaction 
(Sec. IIV Ap and strong dipole-dipole interaction (Sec. 
IIV Bp . In both cases the performance depends critically 
on the Rydberg state lifetime which we calculate for rel- 
evant experimental parameters in Sees. IIV ElIVFl An 
additional aspect of the Rydberg state interactions that 
needs to be addressed is the rate of heating due to differ- 
ences in the ground state and excited state polarizabili- 
ties. We show how to minimize this effect at the expense 
of some additional decoherence in Sec. IjlVDp . 

The results of the calculations of fidelities and decoher- 
ence rates provide a picture of the feasibility of a quan- 
tum logic device capable of executing a large number of 
sequential gate operations. We discuss the expected over- 
all performance of this approach to quantum computing 
in Sec. E] and highlight the areas that are most trou- 
blesome. Possible extensions to the techniques discussed 
here that have the potential for improved performance 
are discussed. 



In its simplest form an optical FORT trap can be cre- 
ated by focusing a single laser beam of wavelength A/ to 
a waist u>/o |33|, |34| • The ground state AC Stark shift due 
to a far-detuned trapping beam is 

U„ = -\\£f\ 2 e* ■ ■ e (1) 

where £f is the amplitude of the optical field, the laser 
intensity is // = eoc|£/| 2 /2, and e is the trapping laser 
polarization. The polarizability is in general the sum of 
scalar, vector, and tensor parts|33 . For a J = 1/2 ground 
state with a linearly polarized trapping beam we need 
only consider the scalar polarizability which we calculate 
numerically using Coulomb wave functions to be ao = 
114 A 3 for a trapping laser at A/ = 1.01 fim. 

The maximum depth of the potential well at the cen- 
ter of the focused beam expressed in temperature units 
is U m — —ao|£/| 2 /4. The spatial dependence of the trap- 
ping potential is then 

e -2( X 2 +y 2 )/™ 2 f (z) 

U f (x,y,z) = U m -3 , (2) 

1 + fr 
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where for a FORT beam propagating along z, wi(z) 



w 



fo 



(l + z 2 /L 2 f), with the Rayleigh length Lf 



/A/. 



For the parameters of Table |U a laser power of 60 mW 
gives U m — 1 mK. 

The FORT can be directly loaded from a MOT pro- 
vided that the product of the capture volume and the 
MOT density is larger than unity. A first approxima- 
tion for the capture volume assumes that all atoms in 
the region where |[///fcs| > T c are captured, while those 
outside this region are lost (in the rest of the paper we 
will express all energies in temperature units and put 
ks = 1)- We expect that the capture temperature T c will 
be similar to T a , the kinetic temperature of the atoms in 
the MOT, provided T a is smaller than \U m \- We define 
a relative trap depth \ = \U m \/T a , so that the capture 
volume vanishes for ^ = 1, and increases monotonically 
with increasing x- A simple calculation then results in an 
expression for the capture volume, 



V 




with z m = Lfy/x — 1- Using numerical values from Table 
Q]we get a capture volume of 8.6 x 10 -9 cm 3 . We have 
found in unpublished experiments that a MOT density of 
a few times 10 8 /cm 3 is sufficient to load single atoms as 
was demonstrated by several groups in recent years [sH 

M 

In the context of quantum logic it is important that the 
atoms are well localized in position and momentum. We 
can estimate the variances of the atomic position and mo- 
mentum by making a parabolic expansion of the FORT 
potential about the origin. The effective spring constants 
of the trap are found to be 



\Um\ 
W% 

\Um\ 

T 2 



(4a) 
(4b) 



and the corresponding oscillation frequencies are 



2 nu„ 

V2 

L f \ in 



1/2 



1/2 



(5a) 
(5b) 



with m the atomic mass. For the above parameters 
and 87 Rb we find w (XjZ) /(27r) = (39,3.6) kHz. At T a = 
50 /iK many vibrational levels will be excited and we 
can use Boltzmann factors to estimate the time averaged 



variances of the position and momentum as 

,,2 



(4) - (vl) - 

(4) = 

(vL) - {<«) - (vL) = 



w 



fo 



4 \u m y 

* 2 w% T a 



2X 2 

Ta 

m 



U„ 



(6a) 
(6b) 
(6c) 



Note that the spatial localization along z can be writ- 
ten in terms of an anisotropy factor = \/ k x / k z = 
TTWfo/Xf such that (z 2 a ) = (x 2 a ) . 

B. Background gas collisions 

Collisions with untrapped background atoms in the 
vacuum chamber result in heating and loss of atoms from 
the FORT and therefore limit the storage time and Tl 
that can be achieved. The characteristic energy change 
for which diffractive collisions must be accounted for is 
~ 2.8 mK for Rbj^]. As we are considering a shallow 
FORT of depth \U m \ = 1 mK we will neglect diffractive 
and heating effects and approximate the FORT lifetime 
due to background collisions mediated by a van der Waals 
interaction as 36] 
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3T b 



in 



-"bCTRb-Rb 



(7) 



with T;, = 300 K the temperature of the thermal back- 
ground atoms of density rib- Using <7R,b-Rb = 2.5 x 
10~ 13 cm 2 |3^| we find t c = 55 s at a pressure of 
1 x 10~ 10 mbar. Even without resorting to cryogenic 
vacuum systems pressures as low as 10~ n mbar are 
achievable, which would imply collisional lifetimes of or- 
der 10 minutes. These estimates are consistent with 
observations [31] of FORT decay times using Cs atoms 
of 50 sec. at pressures of about 10 -10 mbar. 



C. Photon scattering 

Scattering of FORT light by the qubit atoms causes 
some heating and leads to a small amount of decoherence. 
The scattering can be separated into two contributions. 
Elastic or Rayleigh scattering of the FORT light does 
not change or dephase the qubit spin but does heat the 
external degrees of freedom of the atoms. Inelastic or 
Raman scattering occurs at a much reduced rate but, 
since it changes the spin-state of the qubit atom it does 
contribute to decoherence, albeit at a small rate. 

The elastic scattering cross section is 



ki\a Q \ 2 = 1.6 x 10~ 24 



(8) 



where kf — 27r/A/. To get the numerical value we have 
used the parameters given in Table D This scattering 
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produces a heating rate 



dE 



h 2 k 2 f a e I f 



2m fko 



Hu)f<j e If 
2mc 2 



(9) 



Since both the heating rate and the trap depth U scale 
with intensity, their ratio gives the characteristic heating 
time for an atom in the FORT 



\U„ 



dE/dt 



2aom\f 

%CF R 



= 1950 s. 



(10) 



A more conservative definition of the Tl due to Rayleigh 
scattering is the time for the atom to double its motional 
energy, which gives Tl = 97 s for the parameters of Table 

□ 

The inelastic scattering cross section can be expressed 
in terms of the vector polarizability a\ as 



47T , 

T 



= 4^fc 4 K| 2 = 2.3 x 10~ 27 cm 2 



(11) 



where we have used a.\ ~ —6. A 3 . The smallness of the 
inelastic cross section comes from the small coupling of 
the FORT light to the electron spin of the atom which 
scales with the ratio of the fine-structure splitting of the 
Rb P-levels to the detuning of the FORT laser. 

Since the inelastic scattering destroys the qubit state, 
it is a source of decoherence. It is proportional to FORT 
intensity, so it can be reduced if necessary by operating 
at low trap depths. The qubit longitudinal relaxation 
time due to inelastic scattering is 



-Lin — 



1 



3hXj\a \ 



a- m I s /{hid f ) l6w 3 \U m \a 2 
1 



x 151 s, 



\U m \[mK] 

which is very long even for a robust 1 mK trap depth 
D. Laser noise induced heating 



(12) 



laser is Si = 10 12 /Hz. The characteristic time for an 
atom to be heated out of the trap is 



Tl = 



1 



> 2 Si{2i 
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T, 



(14) 



Using the values given in Table IJ and v = 39 kHz gives 
Tl ~ 50 s. If necessary, feedback can be used to reduce 
the laser noise and extend the heating time. 

Fluctuations in the laser beam position are also a 
source of heating. The characteristic heating time can 
be written as 



Tl 



irv 2 S x {v) 



In 



\Urn 



(15) 



where S x is the frequency spectrum of the position fluc- 
tuations. The parameters of Table [I] and Eq. I|6a() 
give yj (aSf) = 0.28 ^m. To obtain Tl = 50 s requires 
\J~S~c = 5.6x 10~ 7 jttm/ \/Hz which is feasible with careful 
attention to mechanical connstruction. 

Finally for an anisotropic trap there is also heating 
from beam-steering fluctuations. This implies, for a 
highly anisotropic trap of aspect ratio £/, a heating time 
of 



Tl 



1 



n 2 v 2 e f SeW. 
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T, 



(16) 



where Sg is the spectrum of angular fluctuations of the 
FORT laser beam. For the parameters we are using £/ ~ 
8 so there is a strong sensitivity to beam-steering noise. 
Nonetheless it should be feasible using fiber optic delivery 
of the trapping beam to achieve Tl ~ 50 s. 

To summarize this section, estimates of storage times 
due to technical noise induced heating are of order 50 s, 
for three different mechanisms. Without appealing to ex- 
traordinary technical developments we can set the total 
contribution due to technical laser noise as Tl ~ 20 s. 
Ultimately this number could be improved by several or- 
ders of magnitude before reaching limits set by quantum 
fluctuations. 



Laser intensity and pointing fluctuations can cause un- 
desirable heating in FORTs[33. The heating rate due to 
intensity noise is 

^± = nVS i (2v)T a (13) 

where v is the trap oscillation frequency, and Si(2v) is 
the one-sided power spectrum of the fractional intensity 
fluctuations. These fluctuations are usually far above the 
shot-noise limit Si = 2hv/P = 4x 10 _18 /Hz at the 1-100 
kHz frequencies of interest here. As indicated above a 1 
mK FORT depth requires 60 mW of At = 1.01 /jm laser 
power which can be readily supplied by a small diode 
laser. A typical fluctuation level for an unstabilized diode 



E. AC Stark shifts 

In the preceding sections we have discussed decoher- 
ence mechanisms that to an excellent approximation af- 
fect the qubit ground states equally. Therefore no de- 
phasing of the qubit basis states is incurred and there 
is no contribution to a finite transverse relaxation time, 
T2. As will be discussed in Sec. IHII we will use the states 
| a) = \F = l,m F = 0) and |6) = \F = 2,m F = 0) as 
our computational basis. In the absence of any applied 
fields these states have a hyperfine splitting U"m(0)//i — 
6.83 GHz. In the presence of a static electric field there 
is a correction to the hyperfine splitting in 87 Rb given 
by|H 
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U hf (E)~U hf (0)\l-a s (0)E 



1 



1 



2(U n= 5 i i = l — J7„ = 5 i l = o) 

I 



U n =i,L=2 — t^i=5,i=0 



(17) 



where a s (0) is the static polarizability and E is the elec- 
tric field amplitude. For a far detuned trapping laser with 
a photon energy that is small compared to the term differ- 
ences that appear in Eq. I|17|l we can estimate the shift by 
making the replacement a s (Q)E 2 — *■ a s (ui f)\£\ 2 /2 so that 
Ea.l(T7)l can be written as Uhf{£) = Uhf(0) — (3U ac where 
/3 = 2Uhf(0) [■■■], with [...] the term in square brackets in 
(H2J. For 87 Rb we find (3 = 4.1 x 1CT 5 . At a trap depth 
of \U m \ = 1 mK the correction to Uu is (3U ac ~ 1500 Hz. 

If the trapping laser had no intensity fluctuations and 
there was no atomic motion this small correction to the 
hyperfine shift would be time independent and would 
only give an unimportant correction to the hyperfine 
splitting. However, intensity noise and atomic motion 
result in a time dependent shift which gives a finite T2. 
We consider first the effect of intensity noise which re- 
sults in state dependent dephasing due to fluctuations in 
Uhf(£)- When the averaged fluctuation vanishes over the 
time scale of interest there will be no additional decoher- 
ence, and indeed intensity fluctuations are not expected 
to be problematic for fis time scale gate operations. How- 
ever as regards storage of quantum information during a 
long calculation it is necessary to consider slow drifts in 
laser intensity that will give qubit dephasing. We define 
an effective T2 due to dephasing by 



T2 = 



2irh 



2ttH I 



6U hf (£) (3u ac sr 



(18) 



where SI /I is the fractional intensity fluctuation. The 
relative intensity noise is a function of frequency. In an 
actively stabilized system the fluctuations will be very 
small at high frequencies. We are most concerned about 
finite fluctuations on the time scale of tens of seconds 
corresponding to the effective Tl given in Table [I] 



The effective T2 is shown in Fig. ^ as a function of the 
fractional laser intensity fluctuation. At a FORT depth of 
| [/„,, | = 1 mK and a relative intensity fluctuation of 10~ 4 , 
which is feasible with active stabilization and well above 
the limit set by quantum noise for mW power FORT 
beams, T2 ~ 12 s. 

Atomic motion within the FORT volume leads to a 
time varying trapping potential, and hence dephasing of 
the qubit states. This problem also arises in the context 
of precision measurements of optically trapped atoms [3^ 
IZoj ]. A rough estimate says that since the atomic po- 
sition spread is approximately (T a /\U m \) 1 ^ 2 Wfo/2, the 
fractional variation in trapping intensity due to the trans- 
verse motion is ~ T a /2\U m \ ~ 0.025 for the parame- 
ters of Table [I] The same fractional variation is also 
found for the axial motion. This implies a motional vari- 
ation in the hyperfine splitting of order 38 Hz. The max- 
imum phase perturbation in one axial vibrational period 
is thus ~ 2n x 38 Hz/3 kHz ~ 0.08 rad. We can also ex- 
press this shift as an effective transverse relaxation time 
T2~ 1/38 Hz = .026 s. 

This time is far shorter than the Tl and T2 times due 
to the other mechanisms discussed above. As has been 
demonstrated experimentally in Ref.|4l| it is possible to 
cancel the differential AC Stark shift of the hyperfine 
states by introduction of a weak beam tuned between 
the hyperfine states that has the same spatial profile as 
the FORT beam. The intensity of the additional compen- 
sation beam can be very low such that the decoherence 
rates due to photon scattering will not change signifi- 
cantly. Kaplan, et al. [4lT | demonstrated a reduction in 
transverse broadening by a factor of 50, and we assume 
that a factor of at least 100 is realistic, in order to arrive 
at the estimate of 2.6 s. given in Tabled Additional dis- 
cussion of motional effects in the context of single qubit 
operations is given in Sec. IIIII 




0.001 0.002 0.003 0.004 

Fractional intensity fluctuation 

FIG. 1: Effective T2 due to differential AC stark shift of 
hyperfine levels for trap depths of 1, 0.5, and 0.1 mK. 



F. Background magnetic and electric fields 

The amount of dephasing caused by trapping and 
background field fluctuations depends on the qubit ba- 
sis states that are used[l3t l39| . As shown in Fig. [2 for 
atoms with nuclear spin 7 = 3/2 there are three possible 
choices of basis states that are first order free of Zeeman 
shifts. In this subsection we consider the three possi- 
ble choices and will conclude that set a) is optimal for 
achieving long storage times with low decoherence. 

The simplest choice (Fig. [3^) is to use \a) = \F = 
l,m F = 0) and \b) = \F = 2,m F = 0). Using the Breit- 
Rabi formula the second order shift of the energy interval 
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FIG. 2: Possible choices of qubit basis states that have a 
magic bias field Bo for which the relative Zeeman shift is a 
quadratic function of field fluctuations. 



expressed as a frequency is 



(gs ~ 9i)^b5B 



1/2 



(19) 



where a>hf is the zero field hyperfine clock frequency be- 
tween the m F = states, /zs is the Bohr magneton, 
gs,gi are the electron spin and nuclear Lande factors, 
and SB is the magnetic field fluctuation. The frequency 
deviation Suiba = ^ba — Whi implies a transverse relax- 
ation time T2 — 2n/8uba which evaluates to T2 = 1740 s 
for SB = 1 mG. While it is in principle possible to shield 
magnetic field fluctuations to an even lower level it is 
difficult to do so in an experiment that requires substan- 
tial optical access to the atom trapping region. Recent 
work|42] has demonstrated suppression of static and fluc- 
tuating magnetic fields to the level of 300 /iG using an 
active feedback scheme. We will assume 1 mG as a con- 
servative estimate of the fluctuation level that can be 
achieved. 

Alternatively we can use the states shown in Fig. 03: 
| a) = \F = l,m F = -1) and |6) = \F = 2,m F = 1). 
At a bias field of Bo ~ 3.23 G the frequency separation 
is quadratically dependent on fluctuations about Bq. A 
fluctuation of 1 mG about the bias point gives T2 = 
2320 s and frequency separations between the qubit basis 
states and neighboring Zeeman states of about 2.3 MHz. 

Finally we could also use the states shown in Fig. [2t: 
| a) = \F = l,m F = -1) and |6) = \F = 2,m F = -1). 
At a bias field of Bq = huj\ 1 {/(2(gs — the fre- 

quency separation is quadratically dependent on fluctu- 
ations about Bq . Defining SB = B — Bq we have 



^>ba 



1 , 4 { (gs-gi)vBSB Y 

3 V ^hf / 



"I 1/2 



(20) 



Apart from a factor of 2/v3 larger sensitivity to fluctua- 
tions we retain the quadratic dependence of Eq. (|19JI at a 
very large bias field. For 87 Rb we find Bq ~ 1219.3 G and 
T2 = 1500 s for SB = 1 mG. At this large bias field the 
separation between neighboring m F levels is hundreds of 
MHz. This large detuning will effectively suppress un- 
wanted transitions during logic operations, but has the 
disadvantage of mixing the hyperfine states so that there 



no longer will be a clean cycling transition between the 
2 S 1/2 \F = 2,m F = 2) and 2 P 3 / 2 \F' = 3,m F > = 3) states 
that can be used for qubit measurement. One possibil- 
ity would be to use a bias field that is applied during 
logic operations, and turned off adiabatically for state 
measurements. 

The use of m F ^ basis states requires that we ac- 
count for m F dependent shifts due to the ground state 
vector polarizability that couples to nonzero ellipticity of 
the trapping laser |39l l43l |44| . The energy shift can be 
written as 



aigpmp 



(21) 



where g F = [F(F + 1) + S(S + 1) - 1(1 + l)]/[F(F + 1)], 
Til p is the total spin projection along the FORT beam 
propagation direction y, and the laser polarization is e = 
(1/V2)(i- V / 1 + e + iyy/l — e). Basis states with mp = 
have no vector shift, and in addition the field insensitive 
states |0) = \F = l,m F = -1) and |1) = |F = 2,m F = 
1) both have gpmp = 1/2 so there is no differential shift 
and no decoherence due to the vector polarizability. On 
the other hand the choice |0) = \F = l,rn F = —1) and 
|1) = \F = 2,mp = —1) leads to a transverse decoher- 
ence time of T2 = (a /ai){I /5I)2-Kh/(k B T m y/i - e 2 ). 
Using the parameters given in the caption of Table 
e = .999 and a fractional intensity fluctuation of SI jl = 
0.001 we find a short coherence time of T2 — 38 ms. We 
see that it is important to avoid decoherence due to the 
vector polarizability so that the preferred choice is the 
qubit states shown in Fig. or [2b . 

Careful analysis along the lines of that used in Sec. 
IIIII shows that the set of Fig. |2b presents significant ob- 
stacles to achieving high fidelity single qubit operations 
using stimulated two- photon Raman transitions. The es- 
sential problem is that the choice of Fig. |2h involves 
ground states separated by |Atof| = 2. For large detun- 
ings the two-photon Raman rate for these transitions is 
proportional to the vector polarizability which has a se- 
lection rule | Ami? | < 1- The Raman rate therefore van- 
ishes so it is necessary to use two Raman beams with 
opposite helicities. In this situation the effective Raman 
rate scales as A e /A^ where A e is the width of the excited 
state hyperfine structure, and Ai is the one-photon de- 
tuning of the Raman beams. However, the Raman beam 
induced ac Stark shifts scale as 1/Ai so that in the limit 
of large detuning the ground state ac Stark shifts be- 
come large compared to the Raman rate between ground 
states. It is not possible, without resorting to more com- 
plex polarization states, to balance the ac Stark shifts of 
the qubit basis states, which leads to entanglement of the 
qubit spin state with the atomic center of mass motion. 
This entanglement represents an undesired decoherence 
mechanism. 

We are thus led to the choice shown in Fig. for the 
qubit basis states. While the mp = states are insen- 
sitive to magnetic fluctuations they are not optimal as 
qubit basis states at low magnetic fields. Atomic motion 
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TABLE II: Physical mechanisms contributing to imperfection 
of single qubit operations. The decoherence probability is 
calculated for a tt rotation, and the fidelity error is calculated 
with respect to an ideal 7r/2 rotation using the metric of Eq. 
ll^f]l . Values listed are for T m = 1 mK, T a = 50 /iK, \ } = 
1.01 /im, Wfo = 2.5 fim, and wo — 5 /im. See text for details. 

in a region of near zero magnetic field is subject to Majo- 
rana transitions between Zeeman sublevels. Transitions 
can be suppressed by Zeeman shifting the states with a 
bias magnetic field. Unfortunately a large bias field con- 
verts the small field quadratic dependence of Eq. H19|) 
to a linear dependence on the field fluctuations about 
the bias point. For example a bias field of 1 G, which 
gives MHz scale Zeeman shifts, with a fluctuation of 1 
mG about the bias point would give T2 = 0.87 s. We 
can do considerably better with a small 15 mG bias field, 
which is sufficient to suppress Majorana transitions, yet 
small enough such that with 1 mG of field fluctuations 
the coherence time is T2 = 56 s. In the remainder of this 
paper we will analyze the implementation of quantum 
logic using the mp — basis states. 

Finally we note that dephasing due to dc electric fields 
is completely negligible. The T2 due to differential ac 
Stark shifts calculated from Eq. I|18[) results from a peak 
electric field at the center of a mK deep optical FORT of 
O(10 6 [V/m]). Since we expect low frequency field fluc- 
tuations to be much less than 1 V/m, the dc Stark shift 
can be neglected. 

III. SINGLE QUBIT OPERATIONS 

A two-site FORT with a single atom loaded in each 
site provides a setting for studying basic one and two 
qubit operations. In this section we start with a study 
of the fidelity and decoherence properties of one qubit 
operations at each site. Of particular concern will be the 
requirement of high fidelity operations at a targeted site 
without unintended disturbance of the neighboring site. 
Simply increasing the separation of the sites to reduce 
crosstalk will imply slow 2-qubit conditional operations, 
so there is inevitably a performance trade-off between 1- 
and 2- qubit gates. We discuss balancing the conflicting 
requirements in Sec. Ivl below. 

Single qubit rotations between ground state levels can 
be performed in several ways. Microwave fields that are 
resonant with Uba can be used, but do not allow direct 



single site addressing. By combining a microwave field 
with an electric or magnetic field gradient a selected site 
can be tuned into resonance The drawback of such 
an approach is that neighboring sites will be subjected 
to off-resonant perturbations. 

Here we analyze an alternative approach using stimu- 
lated two-photon Raman transitions induced by tightly 
focused addressing beams, as shown in Fig. [3| Two- 
photon Raman techniques for laser cooling of neutral 
atoms were pioneered by Kasevich and Chu|46J, and have 
been used recently in optical lattices by the group of 
Jessen0 an d others^]- Raman techniques are also an 
important ingredient in trapped ion experiments As 
the physics of coherent state manipulation with stimu- 
lated Raman pulses is well understood our aim here is 
to analyze a number of contributions to nonideal behav- 
ior that arise in the context of optically trapped atom 
experiments. The physical mechanisms contributing to 
non ideal single qubit operations are summarized in Ta- 
ble [n] and discussed in Sections IIII AlIIl'El In addition 
to qubit rotations fast state measurements are a require- 
ment for error correction in a quantum processor. In 
section All Fl we examine the speed and fidelity of single 
site state measurement using resonance fluorescence. 



A. Speed and decoherence 

As discussed in Sec. [H] the qubit logical basis states 
|a), |6) will be represented with the 87 Rb 5S]/ 2 \F = 
l,mp = 0) and \F = 2,mp = 0) ground state hyperfine 
levels. Qubit initialization to the \b) state can be accom- 
plished by driving the 5S 1/2 (F = 2) - 5P 3/2 (F' = 1) 
transition with a beam linearly polarized along the mag- 
netic field. An additional repumper beam returns atoms 
infrequrntly lost to the F = 1 lower hyperfine manifold. 
When the transition is driven with an intensity several 
times larger than the saturation intensity population will 
accumulate in \F = 2, mp = 0) at the rate of 7/2, where 
7/2tt = 5.98 MHz is the spontaneous decay rate from 
the upper level. The characteristic qubit initialization 
time will be several transfer time constants, or ~ 0.1 fis. 
Experiments with atomic beams |49l| have demonstrated 
preparation purity by optical pumping at the level of 
10- 4 . 

Ground state single qubit manipulations using stim- 
ulated Raman transitions can be performed with high 
fidelity and made quite free of decoherence due to spon- 
taneous emission. The driving fields and atomic level 
structure are shown in Fig. [3J We consider two driving 
fields at frequencies u>\, 102 with detunings Ai = oj\ —u>i a , 
A2 = LO2 — uju, and associated Rabi frequencies fi^a — 
^1,2^1,2/^, with d\ 2 the relevant dipole matrix elements 
between the ground states and the excited state |1). The 
fields propagate along the z axis with polarizations £1,2 
so that the total optical field is 
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FIG. 3: (color online) Single qubit stimulated Raman rota- 
tions. 



where w{z) = w \/l + z 2 / z 2 R , z R = ttw'q/Xr, \ r = Ai, 
and we have taken the Rayleigh lengths of the two fields 
to be equal since |(Ai — A 2 )/(Ai + A 2 )| <C 1. 

Let us assume the atom is in state \a) at t = 0. The 
Raman light is tuned in the vicinity of either the 5P3/2 or 
5P%/2 excited states. For definiteness we consider tuning 
near to the 5Pi/ 2 excited state. In the limit of large single 
photon detuning relative to the excited state hyperfine 
structure the probability for the atom to be in state \b) 
at time t is 



\c b (t)f 



■ sm 



-t 



(23) 



where Q R = Q1CI2/6, 5 = Ai + A 2 , and A = u>\ — u>2 — 
^ a + (|f> 2 | 2 -|^i| 2 )/(2<5). When |A| « \Q R \ « |Ai|, |A 2 | 
we have 8 ~ 2Ai and the effective Rabi frequency is 
£l R ~ niO|/(2Ai). We assume that Ai,A 2 are small 
compared to the 87 Rb fine structure splitting of 7120 
GHz, so we neglect any contribution from the other 5P 
state. 

The probability of spontaneous emission during a tt 
pulse of time t„ = n/Q', with fi' = ^ | O ^ | 2 + A 2 
the effective off-resonance Rabi frequency, is p se — 
t _1 f Q * dt\c p (t)\ 2 , where c p is the amplitude of the ex- 
cited state and r = 27.7 ns is the 5Pi/ 2 lifetime. Ne- 
glecting the two-photon detuning and assuming a piece- 
wise constant pulse profile it is readily shown that p se ~ 
tt/(2|Ai|t). For 87 Rb with Ai/2tt = -100 GHz we get 
Pse = 9 x 10- 5 . 

The speed of two-photon Raman transitions scales 
with the optical intensity. In subsequent sections we 



\F=2>. 




FIG. 4: (color online) Stimulated Raman rotations with two 
excited states. 



will be concerned with corrections to the effective Ra- 
man frequency due to the excited state hyperfine struc- 
ture shown in Fig. 21 We work in the limit of |A| <C 
\Q R \ <C |An|, |Ai 2 |, |A 2 i|, |A 22 |, so the excited states are 
weakly populated. With ei = e 2 = e + the effective Rabi 
frequency is 



|n*| = 



2e 2 (i?5S : 5p)' 
eoch 2 



KI 



A n - A e /4 



2Aii(An - A e ) 



(24) 



with P5S.5P ~ 5.13ao the radial integral between the 
ground and excited states, ao is the Bohr radius, if = 1/9 
is an angular factor, / is the intensity of each Raman 
beam, and A e /2n = 817 MHz is the hyperfine splitting 
of the excited state. Working with 100 /iW in each beam 
focused to a spot with waist wq — 5 /jm at a detuning of 
|Aii|/2tt = 100 GHz we get n R /2n = 4.6 MHz. 



B. Raman beam AC Stark shifts 

In this section we calculate heating and decoherence 
effects due to the tightly focused Raman beams. Before 
proceeding with calculations of the fidelity of qubit oper- 
ations we define the metric to be used. Starting with an 
initial pure state \ip) = c a \a) + Cb\b) a two-photon stim- 
ulated Raman transition results in the transformation 
IV') — > R{Qr, A, t)\ip) where the rotation matrix is 
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TABLE III: Stark shifts of Zeeman ground states due to the Raman beams of Fig. [I] The Stark shifts given by SUr = 
(e 2 R^ S 5Pi ^/(eoch))(K/72)I, with I the intensity of each Raman beam, are evaluated in columns 3 and 4 for / = 100 fiW and 
Au/2tt = -100 GHz. 



zAt/2 



-iAt/2 n R 



n't 

2 



-iAt/2 



ie *A t / 2 r£ sin ^ 



(25) 



In writing Eq. I|25(l we have suppressed a multiplicative 
phase factor and neglected a small correction to some 
of the terms containing A that is proportional to the 
differential ac Stark shift of the basis states. The fidelity 
of a rotation operation i?(0^, A, t) compared to an ideal 
transformation with Rq = R(flnQ, Ao, £o) can be defined 
as 

F = (\(R ^\R^)\ 2 ) , (26) 

where the outer brackets specify an average over any 
stochastic contributions to R. 

In the simplest case of two-photon resonance A = 
and the rotation matrix simplifies to 

' " — [ie 1 ^ sm(6 /2) cos(0/2) J (27J 

where 8 = |fi.jj|i and = arg(ri!/{). For the particular 
case of an ideal ir pulse (6q — it, <fio — 0) with \ip) = \a) 
the fidelity is F = (1/2)(1 — cos^). As this metric is 
independent of errors in the azimuthal angle (f> we find 
it more informative to quantify the fidelity of a single 



In order to give a quantitative account of the ground 
state Stark shifts for Raman beams tuned near the Dl 
line we add the contributions from the \F' = 1,2) excited 
states for both beams to get the results shown in Table 
IIIII Using the physical parameters given in the previ- 
ous section we find the change in trapping potential at 
the center of one of the Raman beams is Ur ~ 220 /iK 
which is a few times less than the \U m \ — 1 mK wells 



I 

operation with respect to an ideal 7r/2 pulse. In this 
case the fidelity is F — (1/2) (1 + cos^sin^). The gate 
errors listed in Table [H] are defined by E = 1 — F. The 
extent to which errors accumulate in concatenated qubit 
operations is an important consideration when designing 
a computational sequence. A discussion of this topic in 
the context of an ion trap experiment has been given in 



Returning to the effect of ac Stark shifts we note that 
in addition to providing controlled rotations between the 
qubit basis states the Raman beams result in unequal ac 
Stark shifts of the Zeeman states which leads to an ad- 
ditional rotation phase 8<p. The rotation phase has an 
average value that must be accounted for[50j, as well 
as a stochastic part resulting from atomic motion that 
leads to a fidelity error. The Raman beam induced Stark 
shifts also play a useful role by enhancing the nondegen- 
eracy of the Zeeman states beyond that provided by the 
very small bias magnetic field. This effect greatly reduces 
leakage out of the computational basis as we discuss in 
Sec. ImDl below. 

I 

created by the FORT beams. Using red detuned Raman 
beams the additional potential is attractive, nonetheless 
the associated dipole forces can lead to heating of the 
trapped atoms. In the limit where the Rabi frequency is 
much larger than the trap oscillation frequency a trapped 
atom will only move a small fraction of its orbit during a 
Raman pulse. The dipole force during the pulse will with 
equal probability accelerate or decelerate the atomic mo- 



10 



tion, so that on average, to lowest order in the ratio of 
trap frequency to Rabi frequency, there will be no heat- 
ing. Another way to see this is to note that the heating 
rate given by Eq. (|13|) vanishes when the perturbation 
has no energy at twice the trap oscillation frequency. 

Nonetheless when we consider the effect of a single 
Rabi pulse there will be a worst case heating of the atomic 
motion that we require to be small compared to the trap 
depth to avoid loss of the trapped atom. It is readily 
shown that the heating due to a single Rabi tt pulse 
is bounded by SU = 2^U R {w 2 j2w 2 ){T a /\U m \){uj x /n R ). 
Using the parameters given in Tables □ and HD and Q R = 
2?r x 4.6 MHz, we find SU ~ 70 nK. Since this amount 
of heating is very small compared to the FORT depth, 
it will not lead to escape of the trapped atom. Note 
that in the limit of fast Rabi frequency the amount of 
heating is proportional to the atomic temperature, but 
is independent of the laser intensity used for the Raman 
pulse. Although in the first approximation the heating 
per Rabi rotation averages to zero there is a contribu- 
tion to the heating rate proportional to the square of the 
trap frequency. For T a = 50 /iK we get an average en- 
ergy increase of 70 nK per operation so the maximum 
number of operations before there is significant heating 
of the atoms is about 10 3 . This implies that recooling 



of the atomic motion after state measurements will be 
necessary to enable many logical operations. 

To quantify the fidelity error due to the Raman beam 
induced ac Stark shifts we note that a Rabi rotation will 
result in a transformation = c a \a) + cj,|6) — * c' a \a) + 
e lS ^c' b \b), where c' a ,c' b are the desired result of the Rabi 
rotation, and S<f> is an additional differential phase shift 
due to the ac Stark shifts. The differential phase can 
be written as 5(f) — Scpb — 8<j> a — t(SUb — SU a )/H with 
t the length of the pulse. The Raman field £ that is 
seen by the trapped atom is time dependent due to the 
atomic motion at finite temperature. The time averaged 
differential phase is proportional to (5(f>) ~ (£ 2 ) and the 
variance of the phase shift is A 2 (S<f>) — (S<f> 2 ) — (Sip) 2 . 
Accounting for the atomic motion in the two transverse 
dimensions, and neglecting the axial motion which gives 
a much smaller contribution to time variation of the field, 
we find (£ 2 ) = £ 2 I [l + (T a /\U m \)(w f0 /w ) 2 ] and (£ 4 ) = 
£q/ [l + 2(T a /\U rn \)(w f /w ) 2 ] , with £ the peak value 
of the field. These expressions are valid in the limit of 
tight confinement where T a <C \U m \. 

Using Eq. and Table UTTI we find to leading order 
in the ratio of the atomic temperature to the trap depth 
for a ir/2 pulse of length t n / 2 = n/2(\il R \) 
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The average phase shift given by Eq. I|28all evaluates to 
0.1 rad for An/27r = —100 GHz. This phase can be com- 
pensated for by adjustment of the relative phase of the 
Raman beams. It turns out using the full dependence of 
differential phase on detuning that there are finite values 
of the detuning such that the ac Stark shifts are equal and 
the differential phase vanishes [&]]]. However, these detun- 
ings are of order the hypcrfinc ground state splitting, and 
are too small to suppress spontaneous emission from the 
5Px/2 states. The standard deviation of the stochastic 
phase given by Eq. (|28b() scales linearly with the factor 
(T a /|C/ m |)(w/o/wo) 2 which expresses the amount of vari- 
ation of the Raman intensity over the cross sectional area 
the atom is confined to. Using Eq. (|2l)|) the gate error is 
E ~ (1/4) A 2 ((50) = 4.4 x 10" 7 at -100 GHz detuning. 



C. Atomic position and velocity fluctuations 

In addition to fluctuations in the phase of the Rabi 
rotation, variations in the atomic position and velocity 
also directly perturb the angle of single qubit rotations 



since the effective pulse area depends on the local value of 
the Raman beam intensities as well as motional detuning 
due to Doppler shifts. Starting with an atom in state \a), 
and applying the Raman fields £\ , £2 the probability for 
the atom to be rotated to state \b) after time t is given 
by Eq. 1)23(1 . Fluctuations in the atomic position and 
momentum lead to fluctuations in the effective pulse area 
at time t n . 

We can characterize a tt/2 pulse and its fluctuations 
due to atomic motion by 

(\QrK /2 ) = 7r/2,(|r! fl | 2 < 2 /2 )=7r 2 /4+(e 2 ),(29) 
= 0, (A 2 4 /2 ) = ( £ 2 ), (30) 

where we have assumed the system has been prepared 
with zero detuning and ei, e 2 are small stochastic param- 
eters. With the Raman beams copropagating the two- 
photon detuning is first order Doppler free. Taking ac- 
count of the velocity spread given by Eq. I|tjc[) we find 
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(4) 
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(31) 



At T a = 50 fiK we find (e 2 .) — 3 x 10~ 13 so we can neglect 
the contribution of Doppler detuning to the rotation error 
and use the simplified rotation matrix of Eq. (|27|l with 
= 0. 

Averaging over the atomic motion in the same way as 
in the previous section we find 



(4) 
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(32) 



and a fidelity error of E = (1/4) (e 2 ) . This result takes 
account of the two-dimensional transverse motion of the 
trapped atom. Adding in the axial motion gives an addi- 
tional factor proportional to £//£ 2 , where £ = ttwo/Xr. 
Our standard system parameters give £//£ 2 = 0.15, so 
this is a small correction which we will neglect. With the 
parameters of Table ITU we find E = 9.4 x 10~ 5 at a tem- 
perature of T a = 50 fiK. This error is larger than that 
due to ac Stark shifts by a factor of (Au/idb a ) 2 , but is 
still very small at typical sub Doppler temperatures that 
are easily reached in a MOT. 

As the motional error scales inversely with the waist of 
the Raman beams it is desirable to use as large a waist 
as possible. In a multiple site device unwanted crosstalk 
occurs if the waist is made comparable to the site to 
site spacing d. Simply increasing d is not feasible since 
that would reduce the fidelity of two-qubit operations, 
as will be discussed in Sec. IIVI The application of Ra- 
man beams giving a Rabi frequency fl R at the addressed 
qubit will result in a leakage Rabi frequency at a neigh- 
boring qubit of fl' R — flne~ 2d ' w ° . For parameters that 
give a 7r/2 rotation at the targeted site there will be a 
fidelity error of E = (tt 2 /16)e~ 4d / w ° at the neigboring 
site. We will use a site spacing of d = 8 /im which gives 
E = 2.2 x 10~ 5 for wq = 5 /im. In practice laser beams 
that pass through a large number of optical elements may 
deviate significantly from an ideal gaussian profile. Full 
characterization of spatial crosstalk will depend on spe- 
cific experimental details. 



D. Polarization effects 

Two-photon stimulated Raman transitions may re- 
sult in the atom being transferred to a Zeeman state 
that lies outside the computational qubit basis. This 
will occur when the Raman beams are a mixture of 
polarization states. The connection between fidelity 
loss due to unwanted transitions and the polarization 
impurity of the Raman beams is calculated in Appendix 
lAl With careful attention to optical design we may 
achieve 8 ~ 10~ 3 for a wide beam. In order to control 
one qubit at a time the Raman beams are focused to 



a waist of wo = 5 /im. Near the focus of a linearly 
polarized Gaussian beam propagating along z the 
positive frequency component of the field can be writ- 
ten as E(x,y,z,t) = e x (£ /2)e- (x2+y2 ^ w2 e^ kz ~ LUt \ 
Consistency with Maxwell's equations requires 
that the actual field is E(x,y,z,t) = [e x — 
i(x/z R )e z ](£ /2)e~^ 2+v2 '>/ w oe< kz ~^ which in- 
cludes a component of e z . Thus a circularly 
polarized field, E ~ e + becomes E(a;, y, z, t) — 
[e+ + {{y - ix)/z R )e ](£ /2)e- ( - x2 +y 2 V w2 e l( - kz -" t l We 
can make a rough estimate of the magnitude of the po- 
larization induced leakage by inserting into Eqs. i|AHA2|) 
\eij\ ~ ^T)/z R = [w f0 /(kw 2 )WT a /\U m \ = 0.0028, 
with the numerical value calculated for our standard 
parameters and T a = 50 /xK. This estimate puts an 
upper limit on the effective polarization purity in the 
interaction region, even when the unfocused Raman 
beams are perfectly polarized. Using this value for 
all coefficients e gives the last column in Table IVII 
which shows that transition amplitudes to undesired 
states will not exceed ~ 6.6 x 10~ 3 . Leakage out of the 
computational basis is a source of decoherence. We 
characterize the decoherence probability for a tt pulse 
by adding the probabilities for leakage out of states \a) 
or | b). We find starting in |o) a leakage probability of 
5.6 x 10~ 5 , and starting in \b) 9.4 x 10~ 5 . We use the 
larger of these numbers in Table HP as an estimate of the 
decoherence probability due to polarization effects. 



E. Laser intensity noise and linewidth 

Intensity fluctuations of the Raman lasers will impact 
the accuracy of the Rabi pulse area. With the aver- 
age intensity and pulse length set to give a ir/2 pulse, 
a relative fluctuation of 5 I/I implies a fidelity error of 
E = (l/4)(<5///) 2 . Active stabilization of the Raman 
laser intensity is limited by shot noise to[l3| 



51 ( Ahjj x 
T " U^/2 



1/2 



where P is the power of the Raman beam and rj is the 
quantum efficiency of the detector in the stabilization 
circuit. With P = 100 fiW, |^|/2tt = 4.6 x 10 6 Hz, and 
T) = 0.5 we find SI/I > 10~ 4 and E > 9.4 x 10~ 8 . 

Finite laser linewidth, and in particular relative phase 
fluctuations of the two Raman beams will lead to errors in 
the phase <p of the qubit rotation. Laser oscillators have 
been demonstrated with a fractional frequency instability 
of 3 x 10~ 16 at 1 sec averaging time [52^. An optical phase 
lock between two laser oscillators with a residual phase 
noise at the 1 /trad level has also been achieved [53. As 
a conservative estimate we will assume the Raman lasers 
can be prepared with a relative phase noise of 1 mrad. 
This implies a fidelity error for a tt/2 rotation of E = 
2.5 x 10~ 7 . 
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F. State detection using resonance fluorescence 

Rapid state selective measurements can be made by il- 
luminating the atom with e + polarized light tuned close 
to the \5S,F = 2,m F = 2) <-> \5P,F = 3,m F = 
3) cycling transition. State readout can be based on 
detection of resonance fluorescence [^3. Alternatively 
amplitude[55j and/or phase shifts imparted to a tightly 
focused probe beam can be used. In either case Pois- 
sonian photon counting statistics result in measurement 
times that are several orders of magnitude longer than 
Rydberg gate operation times. While the use of sub- 
Poissonian light could be advantageous in this context, 
it would add additional complexity. 

It is of interest to estimate the time for performing a 
state measurement with a desired accuracy. The number 
of photon counts recorded in a measurement time r is 



q(r) = r) 



I/Is 



4ir 2 1 + 4 



A 2 



I/I, 



(33) 



where £ld is the solid angle of the collection optics, 7 is the 
radiative linewidth, I is the readout intensity, I s is the 
saturation intensity, and A is the detuning of the readout 
light from the cycling transition. The factor rj < 1 ac- 
counts for the quantum efficiency of the detector as well 
as any optical losses. The counts are assumed poisson dis- 
tributed so that the probability of measuring n counts is 
P n (q) = e~ q q n /nl. We also assume a background count 
rate 60 due mainly to parasitic scattering from optical 
components and detector dark counts that gives a count 
number b(r) = b$T. 

To make a measurement we detect scattered photons 
for a time r. If the number of counts is greater than or 
equal to a cutoff number n c we have measured the qubit 
to be in state \b). The measurement is incorrect if the 
actual number of signal counts was less than n c or the 
number of background counts was greater than or equal 
to n c . The probability of a measurement error is therefore 



E = J2 Pn{q) + £ Pn(b) 
n—0 n—n c 

T(n c ,b) T(n c ,q) 



= 1 - 



r(n c ) F(n c ) 



(34) 



' dt t n - 



where Y(n c ) = (n c — 1)! andr(n c ,g) = / 
the incomplete gamma function. Since T(n c , 0) = T(n c ) 
the error vanishes when b — > and q — > 00. For given 
values of q and b there is an optimum choice of n c that 
minimizes the error. 

Figure [S] shows the error probability as a function of 
measurement time for experimentally realistic parame- 
ters. We see that optimum detection corresponds to a 
very small value for n c and that even with background 
rates as high as b = 10 4 s _1 accurate measurements can 
be made in under 100 /is. A problematic aspect of the 
state measurement process is concomitant heating of the 
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FIG. 5: State detection by resonance fluorescence for two 
different background count rates. The curves are labelled with 
the optimum values of n c . Parameters used were £7^/4^ = .05, 
A/7 = -1/2, I/I s = 1, 7 = 2tt x 6 MHz, and 77 = 0.6. 



atomic motion. The calculations shown assume a detun- 
ing of A = —7/2, so that if counterpropagating readout 
beams are used it should be possible to cool the atomic 
motion while performing the measurement. Experiments 
have demonstrated the feasibility of long measurement 
times exceeding several seconds for single atoms confined 
in micron sized optical traps |32ll5r|. 



IV. TWO-QUBIT PHASE GATE 

In this section we study the performance of the Ryd- 
berg gate using the geometry of two trapped atoms sep- 
arated by a distance R, as shown in Fig. The basic 
idea of the Rydberg gate is to use the strong dipole-dipole 
interaction of highly excited atoms to give a fast condi- 
tional phase shift. The fidelity of a conditional two-qubit 
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FIG. 6: (color online) Geometry of two-qubit interactions. 



13 



operation will be impacted by the mechanisms affecting 
single qubit gates listed in Table[nj plus additional effects 
specific to the use of Rydberg states. As was shown above 
the single qubit imperfections can be quite small, lead- 
ing to projected fidelity errors O(10 -4 ). In the following 
subsections we analyze the additional errors specific to a 
conditional phase gate. 

In comparison to the performance of single qubit oper- 
ations there are two significant complications involved in 
achieving conditional logic. The first is that optical exci- 
tation of Rydberg levels cannot readily be made Dopplcr 
free so that atomic motion introduces pulse area errors. 
The second, and more serious limitation, arises from mo- 
tional heating due to transfer to Rydberg states, and de- 
coherence due to the finite lifetime of the Rydberg states. 
We present a solution to the heating problem based on 
balancing ground and Rydberg state polarizabilities, and 
show that the decoherence rates can be managed such 
that high fidelity gate operations appear possible. 

The Rydberg gate can be optimized in two limits. In 
the first limit (Sec. II V Aj l. where the atoms are rela- 
tively far apart, the two-atom interaction frequency shift 
is small compared to the Rabi frequency of the Rydberg 
state excitation. In the opposite limit of closely spaced 
atoms (Sec. IIVBII the dipole-dipole interaction is large 
compared to the Rabi frequency. In both cases the finite 
lifetime of the Rydberg state sets a lower limit on the 
gate fidelity that scales as the smaller of (VLrt)~ 2 / 3 or 
(Add.T~)~ 2 / 3 , with Qr the Rabi frequency of the Rydberg 
excitation, the dipole-dipole interaction shift, and r 
the excited state lifetime. 

We start the analysis of the phase gate by calculating 
its intrinsic fidelity scaling with speed of operation. The 
dipole-dipole interaction energy is 

1 [ Q 0*i ' 1*12)0*2 - r i2)" 

Mi ' M2 ~ 6 — 



V d d = frAdd = 



4ne Q R 3 



R 2 



(35) 

Here /Xj is the dipole moment of atom j, and 1*12 = Rr 12 
is the atomic separation. When the site to site qubit 
spacing R is sufficiently large the dipole-dipole interac- 
tion shift Add is a small perturbation compared to Qr, 
the Rabi frequency for excitation from \b) — > |r). The 
protocol for a phase gate in this limit is 21]: i) excite 
both atoms from |6) — > |r) with a ir pulse, ii) wait a 
time tdd = 7r/Ad<j(iZ), and iii) transfer both atoms back 
down from \r) — > |6) with a ir pulse. Here Add(R) is the 
dipole-dipole shift at an atomic spacing of R. The result- 
ing idealized logic table is: \aa) — > \aa), \ab) — > — |a6), 
1 6a) — > — 1 6a), |66) — * — 166) which is an entangling phase 
gate. Adding single qubit Hadamards before and after 
the conditional interaction implements a CNOT gate. 
The fidelity of the gate is constrained by the presence 
of four time scales. In the large Rabi frequency limit we 
have 1/uiba <C <C I /Add -C t, where r is the nat- 

ural lifetime of the Rydberg state. The fidelity is funda- 
mentally limited by the combination ojf, a r which should 
be as large as possible. Using characteristic values of 
ujba = 27r x 6835 MHz and r = 100 ^is gives uj^ a T ~ 4 x 10 6 
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TABLE IV: Leading conributions to the fidelity errors of a 
phase gate in the large Rabi frequency limit. See text for 
details. 



which is sufficiently large for high fidelity operation. 

There is an intrinsic source of error in this gate that 
scales with the ratio Add/^R- Assume the Raman fields 
used for the |6) — > \r) transfer provide a ir pulse when 
at most one of the atoms is excited to \r). Then if both 
atoms start in state |6) the excitation to \r) will be im- 
perfect due to the dipole-dipole shift. In this case the 
gate operation will end with a small amplitude for the 
atoms to remain in state \rr). It is readily shown that 
the probability for this to happen, and hence the gate 
error, is E ~ (Add/^R.) 2 ■ This error is intrinsic to the 
design of the gate, so we require a ratio of interaction 
shifts to Rabi frequency of 0.1, or less, for high fidelity 
operation. 

Unfortunately, the finite lifetime of the Rydberg state 
gives a finite decoherence rate, and an error that grows 
linearly with the time of the gate. As the gate time 
scales with I /Add the interaction must be sufficiently 
strong to achieve high fidelity. We can write the prob- 
ability of decoherence as P ~ trfd/r = ir/rAdd where r 
is the excited state lifetime due to all decay mechanisms 
(see Sees. 



IIV El IIV F|) . The total gate error account- 
ing for both imperfect fidelity and decoherence is then 
E ~ (Add/^R) 2 + l/C^Add). The error is minimized for 
A d d ~ (^Ij/t) 1 / 3 , which then gives E - (1/T^ fl ) 2 / 3 . We 
see that excited state lifetime dictates how fast the Ry- 
dberg excitation must be performed to achieve a desired 
error level. 



A. Large Rabi frequency 

In this subsection we investigate operation in the large 
Rabi frequency limit in detail. To quantify the fidelity 
error we average over the four possible initial states of 
the phase gate which gives the results shown in Table 
IIVI The calculations are performed using the rotation 
matrix of Eq. (|25|l and the fidelity definition of Eq. 
(|26|l . The decoherence error listed in the table is the 
integrated probability of a transition out of the Ryd- 
berg state during the gate operation due to spontaneous 
emission or other mechanisms. For example the deco- 
herence error for each atom with the initial condition 
\aa) is E ~ (2/r) J^ 71 "* 1 dt (\n R \ 2 /n' 2 ) sin 2 (ft'i/2) + 
7r/(A dd 2r), with ft' = \/\£Ir\ 2 + uJ 2 a . The error is then 
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doubled to account for two atoms. 

The rotation error corresponds to the probability that 
the atom is in the Rydberg state at the end of the gate. 
For the initial condition \bb) the two atom error is 



E = 1 



J) R(n R ,& dd ,K/\n R \) (J \ 



R{n R> A dd ,n/\fl R \) 



(36) 



Here the rotation matrices operate on the states |6) and 

|r) and the vector (l 0) represents \b). The last row in 
Table llVl gives the gate error averaged over the possible 
input states. While a particular computation may weight 
certain input states more heavily than others the average 
gate error is indicative of the gate performance. The 
average error is minimized for 



*dd I opt 



4r 



1/3 



4w, 2 

ba 



1/3 



which determines the leading order contributions to the 



error as 



EL 



Ifi/ 



■pt 



3 
2V3 



rin, 



2/3 



(37) 



The optimum dipole-dipole shift, gate time gate time 
i g = 2ir/\Q R \ + 7r/ A dd | opt, and gate error are shown in 
Fig. [7| We see that errors less than 10 -2 are possi- 
ble with Rabi frequencies of a few tens of MHz and a 
gate time of less than 1 /is. The corresponding optimum 
dipole-dipole shifts of order 1 MHz can be achieved at 
tens of microns of separation as we discuss in Sec. IIVCI 

Finally we note that the absolute minimum gate er- 
ror occurs for \£l R \ intermediate between u>b a and 1/r. 
Minimizing Eq. I|37|) with respect to we find 

l^lopt = [ 7rw & a /(2 1 ^ 2 ' r )] 1 ^ 4 and a minimum error of 
-E|min = [2 7 / 2 7r/(Tcjb a )] 1 / 2 . For the parameters used in 
Fig. Qwe get 1^1 opt = 2?r x 183 MHz and E\ min = 
2.9 x 10" 3 . 



1. Additional imperfections 

There are two additional imperfections when working 
in the large Rabi frequency limit. The first of these is 
due to fluctuations in the distance between the atoms. 
For atomic separation large compared to the extent of 
the thermal motion the actual qubit separation will be 
R + Xl +x 2 , with (x 2 ) = {xf) = {w 2 /A){T a /\U m \). The 
variance of the interaction phase is easily shown to be 
A 2 ((f>) = {9n 2 /2)(w 2 f0 /R 2 )(T a /\U m \). Using this variance 
as a measure of the gate error gives E = (l/4)A 2 (</>) x 
(1/4) = 3.5 x 1(T 4 at R = 50 (j,m. The final factor of 1/4 
accounts for an averaging over the possible initial states. 
We see that at R — 50 fim this error is small compared to 
the intrinsic gate error which dominates for R > 15 /zm. 
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FIG. 7: (color online) Two-qubit phase gate performance in 
the limit of large Rabi frequency: a) optimum dipole-dipole 
shift and gate time, and b) minimum achievable fidelity error. 
Calculations for r = 100 tis and ujba =2ix 6835 MHz. 



The second subsidiary imperfection is the presence of 
heating due to two-body forces when both atoms are 
excited to the Rydberg state. This effect could limit 
the number of gate operations before motional cooling 
is needed and lead to decoherence through undesired en- 
tanglement of the motional and spin states. The heating 
rate can be estimated simply as P — Fv ~ Ftofuufo. 
Using F = —dVdd/dR for the two body force we get 
a peak heating power of P ~ 3hAddUJ x w /o / R- Using 
Add = 2n x 1 MHz and R = 50 /xm gives P ~ 1.8 fj,K//j,s 
for our standard FORT parameters. We therefore expect 
a maximum of about 1 fiK of heating when both atoms 
are initially in the state \b) which is coupled to the Ry- 
dberg states. Although the heating power will average 
to zero over many operations there is a finite probabil- 
ity for undesired motional entanglement. The spacing 
of the radial vibrational levels in temperature units is 
Ai^vib = hu) x /kB — 1.9 /iK which is comparable to the 
peak heating value. On average there will be a reduced 
probability for a change of the vibrational state since the 
atoms spend proportionately more time near the turning 
points of the motion where the velocity is small. 

The above errors due to fluctuations in the atomic sep- 
aration and two-body forces are specific to gate opera- 
tion in the large Rabi frequency limit. There is an ad- 
ditional error source that is common to both modes of 
gate operation, which is the presence of Doppler shifts 
of the Rydberg excitation beams due to atomic motion. 
For ground state qubit rotations two-photon stimulated 
Raman transitions are essentially Doppler free for co- 
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TABLE V: Leading conributions to the fidelity errors of a 
phase gate in the large dipole-dipole shift limit. See text for 
details. 



propagating beams, and the error given by Eq. I|31|) is 
insignificant. Two-photon excitation of Rydberg levels 
using 0.78 fim and 0.48 /urn beams as indicated in Fig. 
cannot be made Doppler free. For T a = 50 fj,K and 
\Q R \ = 2tt x 10 MHz we find using Eq. J3IJ, fidelity 
errors of 3.3 x 10 -4 and 1.9 x 10~ 5 for co- and counter- 
propagating beams respectively. These errors are signif- 
icantly smaller than the intrinsic gate errors shown in 

Figs. EE 



B. Large dipole-dipole frequency shift 

There is an alternative mode of operation of the phase 
gate that removes the dependence on variations in the 
interatomic separation and eliminates the two body heat- 
ing discussed above. This mode can be used for closely 
spaced atoms for which 3> \Qr\- In this limit 

the time scales have the ordering l/wba <C I /Add <C 
l/|f2iil <C t. The protocol for a conditional phase gate is 
then |2lJ: 1) excite atom 1 (control atom) from \b) — > \r) 
with a 7T pulse, 2) excite atom 2 (target atom) from 
1 6) — > \r) — > \b) with a 2-7T pulse, and 3) transfer atom 
1 back down from \r) — ► \b) with a ir pulse. Note that in 
contrast to the protocol used in the large Rabi frqeuency 
limit we assume here that the atoms can be individually 
addressed. 

The resulting idealized logic table is: \aa) — > \aa), 
\ab) — > \ab), \ba) — > \ba), \bb) — > — \bb) which is an en- 
tangling phase gate. The gate errors are calculated using 
the same techniques as in Sec. IIV Al which results in the 
errors shown in Table |V| The average error is minimized 
for 



ill opt 



1/3 



which determines the leading order contributions to the 
error as 



E 



3tt 2 / 3 



1 



opt 



25/3 ( Add r)2/3 



7A dd 

3-L 



(38) 



The optimum Rabi frequency, gate time t g = 
47r/|Oij| opt |, and gate error are shown in Fig. |S] We 
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FIG. 8: (color online) Two-qubit phase gate performance in 
the limit of large dipole-dipole shift: a) optimum Rabi fre- 
quency and gate time, and b) minimum achievable fidelity 
error. Calculations for r = 100 fj,s and Uba ~ 2tt x 6835 MHz. 



see that errors much less than 10~ 2 are possible with 
dipole-dipole shifts of a few tens of MHz and a gate time 
of less than 1 /is. This protocol appears more promising 
for achieving very small gate errors than the large Rabi 
frequency protocol, since it is easier to achieve very large 
dipole-dipole shifts than it is to achieve very large Rabi 
frequencies. 

Finally we note that minimizing Eq. I|38|) with respect 
to Add we find Add\ op t = yS/liuba and an absolute 
minimum error of E\ min = [170l7r 2 /(128r 2 w 2 Q )] 1 / 3 . For 
the parameters used in Fig. ||we get A^opt = 2n x 
3160 MHz and E\ min = 1.9 x 10~ 4 . 



C. Dipole-dipole interaction strength 

The long-range interactions between two Rydberg 
atoms are extremely strong and are the heart of the quan- 
tum computation scheme discussed in this paper. We 
consider here the Rydberg-Rydberg interactions in two 
limits: zero electric field, where the long-range 1/R 6 van- 
der-Waals interactions are of primary importance, and 
in a hybridizing electric field, where the atoms attain a 
"permanent" dipole moment of magnitude n 2 eaQ. 

For S states, the van-der-Waals interaction is domi- 
nated by the near resonance between the states nS + nS 
and nP+ (n — 1)P, as can be seen from Fig. [5] In the ap- 
proximation that we neglect the contribution from other 
states, the Rydberg-Rydberg potential energy curves are 
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simply given by 



V(R) = - 2 



4W f 
3 4 



(39) 



where U 3 (R) = e 2 (50S\\r\\50P}(50S\\r\\49P}/(4ire R 3 ) = 
5.75 x 10 3 MHz ^m 3 /R 3 and the SS-PP energy defect is 
S = E(49P) + E(50P) - 2E(50S) = -3000 MHz. At the 
10 /im separations of interest here, U3 <C 5 so that 



V(R) 



35 



(40) 



As can be seen from Fig. the van-der-Waals potential 
for 505* is insufficient to allow for the MHz processing 
that is desired here. 

Van-der-Waals interactions between D-states are dom- 
inated by the near-resonance of nD + nD and (n + l)P + 
(n — 1)F. In general, the van-der-Waals interactions are 
of comparable strength when compared to the S-states. 
Because of the degeneracy of the D levels, however, it 
can be shown that U3 vanishes for one each of the 3 S+ 
and 1 S+ molecular states [57|- Thus there is no blockade 
at zero field for these states. 

The scaling of the van-der-Waals interaction with n 
is very rapid, V oc n 11 . For n — 100, V(10 /im) ~ 50 
MHz. While this is sufficient for a conditional phase gate 
(see Fig. [8J we can considerably enhance the interaction 
strength by applying a dc field. Thus we now consider 
the long-range interactions in the presence of an applied 
electric field. 

There are two primary ways to enhance the Rydberg- 
Rydberg interaction using an applied electric field. The 
first is to use the electric field dependence of the state en- 
ergies to bring the energy defect 5 to zero, in which case 
V{R) ~ 2U?,{R) I '\/3. This would produce an extremely 
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FIG. 9: (color online) Stark map for states near n = 50 for 
87 Rb. At electric fields around 5-6 V/cm, the atom acquires a 
large permanent dipole moment, oriented in space along the 
applied field. 
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FIG. 10: (color online) Isotropic dipole-dipole interaction for 
excitation to the Rb 50S state. The energy defect between 
the 50S+50S and 49P+50P states significantly reduces the 
interaction as compared to the ideal degenerate case. Going to 
larger n partially compensates for the non-zero energy defect. 



strong isotropic interaction of maximum strength, for ex- 
ample n — 95 gives a frequency shift at 10 /jm separation 
of 160 MHz. However, inspection of Fig.[5]shows that the 
P-states tune the wrong way in an electric field, increas- 
ing rather than decreasing 5. For D-states there remains 
the problem of zero dipole-dipole interactions for some 
molecular symmetries. 

The second way to adjust the Rydberg-Rydberg in- 
teractions is to use a large enough field to strongly mix 
states of different L. In this case the atom acquires a 
field-independent dipole moment, as suggested for exam- 
ple by the linear field dependence of the 50D state at 
electric fields between 5 and 6 V/cm. In this case the 
Rydberg-Rydberg interaction is strong but anisotropic 
and we have the dipole-dipole interaction of Eq. (|35|l 



V dd (R) = HA dd (R) 



4ire R 3 



[l-'Scos 2 9) (41) 



where 9 is the angle between the interatomic axis and the 
electric field. 

In the proposed geometry for the quantum computer, 
the electric fields can be aligned to 9 — for maximum 
interaction strength. The electric dipole moment is on 
the order of n 2 eao, giving a dipole-dipole interaction of 
comparable size as the zero field case with 6 = 0. We have 
numerically estimated the dipole moments for the field- 
mixed D-states, and obtain for example fi = 3300 eao for 
n = 50. The resulting interaction strengths as a function 
of distance are shown in Fig. El We see that interac- 
tion frequencies in excess of 100 MHz can be achieved 
for R = 8 /im at n = 70 and well beyond R = 10 /j,m 
at n = 95. Thus the dipole-dipole interaction strength 
needed to optimize the phase gate, as shown in Fig. [SI 
can be achieved for qubit separations tht are optically re- 
solvable. It has been pointed out|53 that n should not be 
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FIG. 11: (color online) Dipole-dipole interaction strengths as 
a function of interatomic distance for various principal quan- 
tum numbers. 



too large in order to avoid collisions between the valence 
electrons of neighboring Rydberg atoms. For R = 8 /im 
this implies n should be less than about 100. Although 
n = 100 is sufficient for the geometry considered here, 
it is also true that in the large dipole-dipole shift limit 
discussed above this limitation is not needed as only one 
atom at a time is excited to a Rydberg level. 



D. Balancing the ground and excited state 
polarizabilities 

The implementation of the two qubit conditional gate 
described so far has an inherent flaw that stems from 
the need to make real transitions to Rydberg levels. The 
atoms are confined in attractive optical potential wells 
since the ground state polarizability is positive for light 
tuned to the red of the first Rb resonance lines. How- 
ever, for the same red detuned trapping laser, the high 
lying Rydberg levels have a negative polarizability that 
provides a repulsive potential. Excitation of an atom 
to the repulsive Rydberg state during a gate cycle leads 
to heating and decoherence through entanglement of the 
spin and motional states. Looking at Fig. ^| we see 
that for Xf = 1.06 /zm the 50_D 5 / 2 Rydberg level has 
a polarizability of ao(50D 5 / 2 ) = —77 A 3 . This polar- 
izability is about 95% of the free electron polarizability 
of —e 2 /m e uj'j, with m e the electron mass. This implies 
that the Rydberg polarizability is only weakly dependent 
on the principal quantum number n for large n, so the 
discussion given here for the case of the 50Z? 5 / 2 level is 
applicable to any of the highly excited states. 

When the Rydberg state sees a repulsive potential the 
heating is minimized by turning off the trapping laser 
during the gate operation. For atomic temperatures large 
compared to the trap vibrational energy the motion is 
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FIG. 12: (color online) Polarizability of the 5 Si/ 2 and 50D 5 /2 
states. The inset shows the behavior near the 5S — 6P tran- 
sition. 



quasi-classical and the average heating per gate cycle of 
duration t is (At/) = T a ^(ujT) 2 , where w is the radial 
trap vibration frequency. For r = 0.5 /is and uj = 2tt x 
39 kHz we get (AU) /T a = 7.5 x 10~ 3 . This heating rate is 
about 5 times larger than that due to the AC Stark shifts 
of the Raman lasers discussed in Sec. 



IIII Bl and implies 
that several hundred operations can be performed before 
there is significant heating of the atomic motion. 

The heating can be eliminated by choosing a FORT 
laser wavelength that gives equal polarizability for the 
ground and Rydberg states. Referring to Fig. one 
possibility is to tune the FORT laser to the blue side of 
the BOD — 6P transition where the Rydberg level acquires 
a positive polarizability. The scalar polarizability of the 
Rydberg level is given by 



a (nD 5/2 ) 



1 



J 7' J',-yJ 



3h2J 



1 j/ji U *>J'„J 



-\(l'J'\\D\\lJ)\' 
(42) 

where D is the electric dipole operator and fouj^i ji ~j = 
Ejiji — E-yj. Away from the resonance the polarizability 
changes slowly with ujf and for the |50D 5 / 2 ) level it is 
ao(50_D 3 / 2 ) — —69 A 3 . When ojf is near resonant with 
the \nD 5 / 2 ) 1 6-P3/2) transition the polarizability is the 
sum of the background value ao plus the resonant con- 
tribution to the sum. Using Coulomb wavefunctions for 
n = 50 we find (50Z5 5 /2 1 1 ^ 1 1 / 2 ) = 0.059 ea so that a 
rough estimate for the detuning condition at which the 
polarizabilities are equal is 



A f = 



1 1 



|(50£> 5/2 ||£>||6P 3 /2> 



3h 2 J + 1 a (5S 1/2 ) 
= 2ir x 945 MHz. 



00(50,0, 



5/2) 



(43) 



Unfortunately there is a penalty associated with bal- 
ancing the polarizabilities in this way since there is a 
probability for the gate operation to end with the atom 
having finite amplitude to be in the I6P3/2) level which 
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lies outside the computational basis. The Rabi fre- 
quency for this transition (starting from mj = 3/2) 
is given by = (£ f /(hVl5))(50D 5/2 \\D\\6P 3/2 ) and 
the probability for population transfer to the |6P 3 / 2 ) 
state at the end of a Rydberg operation is bounded by 
P max = \Q f \ 2 /(\fl f \ 2 + Aj). For a trap depth of 1 mK 
we find |%| = 2vr x 210 MHz and P max = 0.16. This 
upper bound assumes that the product of the gate time 
and the spontaneous lifetime of the 6P level is not large, 
otherwise there is an additional loss mechanism due to 
decay out of the 6P level. 

This decoherence probability is roughly proportional to 
the FORT laser intensity and can be reduced by working 
with a shallower FORT trap. This suggests that atoms be 
loaded into 1 mK deep traps where they can be cooled to 
sub Doppler temperatures with standard methods, fol- 
lowed by adiabatic reduction of the trap depth before 
performing logical operations. With T a = 5 /jK and 
\U m \ = 100 fiK the atom localization and motional errors 
will be the same as considered in the rest of this paper, 
while the decoherence probability per Rydberg operation 
will be bounded by P max = 0.016. 

While we expect this to be a viable approach for initial 
experiments it can be readily shown that the decoherence 
probability with polarizability balancing of a nD Ryd- 
berg level scales proportional to 1/| {^^5/2! I-Dl I6P3/2) | 2 - 
Since |(nP 5 /2||-D||6P3/2)| 2 ~ l/n 3 the leakage problem 
scales as n 3 and the technique is most useful for low lying 
Rydberg states that have large transition dipole moments 
with the 6P level. 

An alternative solution is to choose a FORT wave- 
length such that the ground state acquires a negative 
polarizability which coincides with that of the Rydberg 
level. In this situation the optical potential is negative so 
the atom is trapped at a local minimum of the intensity. 
This approach has a number of advantages since for an 
intensity profile of the form / = I (l - e~ 2r the 
time averaged trapping intensity at the position of the 
atom is T a /\U m \ times smaller than it would be for an 
attractive potential of the same depth. This increases the 
Tl and T2 times in rows 2,3 and 6 of Table [I] associated 
with photon scattering and motional AC Stark shifts by 
the same factor. It was pointed out by Safronova et al. 
|59|| that the polarizabilities can be balanced in this way 
by tuning the trapping laser between the D lines. For 
87 Rb the balancing point is at \f ~ .79 fim. However 
the vector polarizability at this point has a very large 
value of a>i ~ 3600 A 3 which would lead to a Tl due to 
inelastic scattering (see Eq. 1121 of about 1 ms for the 
parameters of Tabled which is unacceptably short. 

A more favorable alternative is to use a FORT wave- 
length of A/ ~ 0.431 /im which gives equal polarizabilities 
of ao ~ —12.5 A 3 , as shown in the inset to Fig. At 
this wavelength a± ~ —0.01 A 3 which gives an inelastic 
scattering Tl that is longer than any other time scale in 
Table [i] The only drawbacks of this solution are techni- 
cal, not fundamental. There is additional experimental 



complexity associated with creating optical beams with 
local intensity minima, as well as the need for a medium 
power laser in the blue. Recent developments in solid 
state laser sources and parametric frequency convertors 
render the latter requirement readily solvable. 



E. Rydberg state radiative lifetimes 

The lifetimes of highly excited Rydberg nl states are 
affected strongly by background blackbody radiation. If 

the K lifetime is T^r the finite temperature lifetime can 
be written as 



1 
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(44) 



where 



r (bb) 



is the finite temperature blackbody contri- 



bution. The K radiative lifetime can be calculated by 
summing over transition rates or approximated by the 
expression 60J 
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(45) 



For all the alkalis ~ 3. Parameters for Rb are as 



2.94, a P = 3.02, a D = 2.85, a F = 2.95, and r 



(o) _ 



1.43, 



r p 0) = 2.76, = 2.09, and - 

For large n the blackbody rate can be written approx- 
imately as 61] 
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where a is the fine structure constant. Equation (|46ll in- 
cludes transitions to continuum states so that it accounts 
for blackbody induced photoionization. Figure IT3l shows 
the radiative lifetime t„l for n up to 100 and several I 
states. We see that for n > 65 the S,P,D, and F states 
have lifetimes greater than 0.1 ms at room temperature. 



F. FORT trap induced photoionization 

Highly excited states are also unstable against pho- 
toionization from the intense trapping light. Since the 
Rydberg electron is nearly free, and the photoionization 
is far above threshold, the photoionization cross sections 
are small. We have estimated the photoionization rate 
of high Rydberg states by calculating the cross section 
using Rydberg and continuum wavefunctions. 

The cross section is given by [(H) 



2tt 



, He 2 df_ 
mc dE 



(47) 



where / is the oscillator strength and E c = Hlu + E r the 
energy of the continuum electron produced by the photon 
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FIG. 13: (color online) Excited state lifetime due to ra- 
diative decay for T = (blue, thin lines) and T = 300 K 
(red, thick lines) for S,P,D, and F states of Rb. 



of frequency u> from the Rydberg state of energy E r . The 
oscillator strength distribution is 



dE ^ 



M(2L r 



1 



^ r (r)r(p LtE (r)dr 



(48) 



where L> is the larger of L r and L and the continuum 
wavefunction is normalized per unit energy 



4>l,e{ t ) 



2m 



sin(kr + S l ,e)- 



(49) 



We have used two methods to detemine the bound and 
continuum wavcfunctions. In the first we use quan- 
tum defect theory^ to find the Coulomb wavefunctions 
and phase shifts S^e- The second method uses the In- 
dependent model potentials of Marinescu et al. |r33| , with 
values slightly adjusted to give the proper quantum de- 
fects of the Rydberg states. The Numerov method was 
then used to find 4>l,e- Numerical results for the two 
approaches agree within a factor of about 2 for the range 
of n discussed here. The discussion and Fig. ^]give the 
results obtained using the Marinescu potential. 

The n = 50 photoionization cross sections are shown as 
a function of L in Figure [21 The photoionization cross 
section for the S-state is much smaller than L > 1 states 
due to the it/ 2 phase shift between the S and P wave- 
functions pl|. The cross sections increase dramatically 
for the P and D states before slowly decreasing with fur- 
ther increases in L, While the 50S cross sections are very 
small, the 100 x larger cross section for the higher I lev- 
els implies that the photoionization rate for the S— states 
will depend sensitively on mixing with the P— levels due 
to external fields. Thus the photoionization cross sec- 
tions may be as large as 10 -20 cm 2 , or as high as 7 P i = 
31,000/s for a 1 mK trap depth. 

The photoionization cross sections decrease with n, a 
factor of 6 going from n = 50 to n = 90 as shown in 
the inset to Figure El Thus at n = 90 we get a mini- 
mum lifetime of r p i ~ 190 fis for the P and D states and 
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FIG. 14: (color online) Photoionization cross section vs 
ejected electron energy, for the Rb 50S Rydberg state with 
Xf = 1.01 /Ctm. 



substantially longer for the other L states. Taking into 
account the blackbody lifetime calculated in the previ- 
ous section we conclude that the Rydberg lifetime will 
exceed 100 [is for n > 80. This validates the gate fidelity 
estimates discussed in Sees. IIV Al and IIVBI Since the 
photoionization rate scales linearly with the trap depth 
even better performance is possible by further cooling of 
the atomic motion and a corresponding reduction of the 
trap depth, or by use of a blue detuned FORT laser as 
discussed in Sec. IIVDI 



V. DISCUSSION 

In this paper we have presented a detailed analysis of 
quantum logic using neutral atoms localized in optical 
traps. Two-photon Raman transitions are used for single 
qubit gates, dipole-dipole interactions of Rydberg states 
provide a two-qubit conditional phase gate, and detection 
of resonance fluorescence is used for state measurement. 
We have concentrated on an implementation of neutral 
atom gates that should be feasible with currently avail- 
able experimental methods and laser sources. The Tl 
and T2 coherence times for qubit storage, and error es- 
timates for single qubit gates are summarized in Tables 
[3 and [H] The intrinsic two-qubit gate errors are shown 
in Figs. and [SI The analysis supports the feasibility of 
MHz rate logical operations with intrinsic errors O(10 -4 ) 
for single qubit operations, O(10~ 3 ) for two qubit gates, 
and state measurements in less than 100 /is with O(10 -3 ) 
measurement error. We show that coherence times of at 
least several seconds are possible in red-detuned attrac- 
tive optical traps. Taken together these numbers sug- 
gest an attractive framework for experimental studies of 
quantum logic with neutral atoms. 

It should be emphasized that the Rydberg gate ap- 
proach does not rely on cooling the atoms to the mo- 
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tional ground state. While ground state cooling has been 
demonstrated in optical traps, and lower atomic temper- 
atures will lengthen coherence times and reduce some of 
the gate errors, the difficulty of maintaining the atoms 
in the motional ground state in the presence of heating 
mechanisms should not be overlooked. Our assumption 
of T a = 50 ^iK does not require complex cooling schemes, 
and implies that O(10 3 ) single qubit logical operations 
can be performed without significant reduction in fidelity 
due to motional heating. 

As discussed in Sec. IIV Dl the differential polarizabil- 
ity of the ground and Rydberg states in a red detuned 
FORT will lead to substantial heating or loss of coher- 
ence. While initial experiments are viable in red detuned 
FORTs our analysis suggests that a blue FORT where the 
atoms are trapped at a local minimum of the intensity 
will ultimately be necessary to realize the full potential 
of this scheme. The blue FORT will also substantially 
improve the coherence times for qubit storage, and to a 
lesser extent the Rydberg state lifetime. 

Extending the two-qubit approach described here to a 
large number of qubits will involve solving challenges re- 
lated to loading and addressing of multiple sites. The 
Rydberg gate approach does appear intrinsically well 
suited for implementation in a two-dimensional array, in- 
cluding error correction blocks. The large dipole-dipole 
shift limit can be used for gates between neighboring 
sites, while the large Rabi frequency limit which works 



at longer range may allow non-nearest neighbors to be 
coupled. Figure [7| shows that in the limit of large Rabi 
frqeuency, gate errors less than 10 -2 are possible with 
dipole-dipole coupling strengths of only 1 MHz. For 
Rydberg levels with n = 70 we can achieve a coupling 
strength of f MHz at a separation of about 40 /im. With 
a 8 /im qubit spacing this suggests the possibility of cou- 
pling blocks of 25 or more qubits without physical mo- 
tion. By taking advantage of the directional properties of 
the dipole-dipole interaction described by Eq. (|35|1 it is 
possible to perform row parallel operations between pairs 
of qubits with strongly suppressed crosstalk. While there 
are many appealing features of the approach studied here 
we emphasize that the extent to which it will prove pos- 
sible to perform arbitrarily large, scalable quantum com- 
putations, remains an open question that will require a 
great deal of further theoretical and experimental stud- 
ies. 
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APPENDIX A: RAMAN TRANSITIONS DUE TO 
IMPERFECT POLARIZATION 



We account for imperfect polarization of the Raman 



beams by putting e\ = (e_| 



ei_e_ 



■eioe )/(l 



elo) 1 / 2 and e 2 = (e+ + e 2 -e_ + e 20 e )/(l + e 2 2 _ + e^) 1 / 2 . 
Referring to Fig. 0] and limiting ourselves to transi- 
tion amplitudes that are linear in the small parameters e 
we must account for the following undesired transitions 
starting from \a) = |10): 

| a) — > 1 1 — 1> eio(£ie )(£ie + )* + e 20 (£ 2 e )(£ 2 e + )* 

\ a ) -> I 11 ) e *o(^i e +)(^i e o)* + e 2 o(£2e+)(£ 2 e )* 

|o> — ► |2 - 2) ei _(£ie_)(£ 2 e + )* 

|a)-»|2-l) ei (fie )(£ 2 e + )* 



|a) -> |21) e; o (£ ie+ )(£: 2 eo)* 
|a) -> |22) e|_(£ ie+ )(£ 2 e_r 



(Al) 



We have labelled the kets as \F, mjr) and the factors on 
the right indicate the participating fields and polariza- 
tions. The Clebsch-Gordan factors and detunings associ- 
ated with these transitions are different than those that 
apply for the desired \a) — > \b) transition. We will as- 
sume that the fields are in two-photon resonance for the 
| a) — > | b) transition accounting for the ac Stark shifts 
given in Table ITTTI The Rabi frequencies normalized to 
that of the desired transition are given in Table IVT1 which 
shows the transition amplitude for excitation of an un- 
desired state with parameters that give a n pulse on the 
desired transition. 
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2.42 
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-2.41 
2.42 
-4.79 
-2.39 
2.38 
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TABLE VI: Transition amplitudes out of the computational basis due to imperfect optical polarization. The effective Raman 
frequencies are given by Qr = (2e 2 B2,s,5P 1/2 I (toch 2 ))(K/72)I , with I the intensity of each Raman beam. The last column 
shows the amplitude of the leakage state for a n pulse on the \a) — > \b) (or \b) — > \a)) transition, evaluated using Eq. I23H with 
/ = 100 /jW and Aii/2tt = -100 GHz such that \n R (\a) -> \b))\/2n = 4.6 MHz. 



